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ABSTRACT 

The purpose of this investigation is to develop a new model for the acceleration 
of the fast solar wind by nonlinear, time-dependent multidimensional MHD 
simulations of waves in solar coronal holes. Computational studies indicate 
that nonlinear waves are generated in coronal holes by torsional Alfven waves. 
These waves in addition to thermal conduction may contribute considerably to 
the accelerate the solar wind. Specific goals of this proposal are to investigate 
the generation of nonlinear solitary-like waves and their effect, on solar wind 
acceleration by numerical 2.5D MHD simulation of coronal holes with a broad 
range of plasma and wave parameters: to study the effect of random disturbances 
at the base of a solar coronal hole on the fast solar wind acceleration with a 
more advanced 2.5D MHD model and to compare the results with the available 
observations; to extend the study to a full 3D MHD simulation of fast solar wind 
acceleration with a more realistic model of a coronal hole and solar boundary 
conditions. The ultimate goal of the three year study is to model the fast solar 
wind in a coronal hole, based on realistic boundary conditions in a coronal hole 
near the Sun. and the coronal hole structure (i.e., density, temperature, and 
magnetic field geometry) that will become available from the recently launched 
SOHO spacecraft. We report on our study of the effects of slow magnetosonic 
waves recently detected by EIT on the solar wind. 


1. Introduction 

The solar corona and the solar wind play a major role in determining the physical 
conditions of the interplanetary environment. However, the exact mechanism that energizes 
the solar coronal plasma and the solar wind is poorly understood. To study the physics 
of the multi-species plasma in the inner corona (the region where most of the heating and 
acceleration takes place) we use the multi-fluid description, applicable dose to the sun. 
This is the first multi-fluid, multi-dimensional study of the solar wind that includes the 
self-consistent effects of MHD waves in a model coronal hole near the sun. 

Recent Solar and Heliospheric (SOHO) Ultraviolet Coronagraph Spectrometer (UVCS) 
observations show that the physical properties of the minor ions in the corona provide clues 


on the coronal heating and solar wind acceleration mechanism. UVCS observations suggest 
that protons roach 3.5 x l O 1 * K within several solar radii, and 0~ :i and other minor ions 
exceed temperatures of 10' K in coronal holes (Kohl et. al. 1997). 

Low and high frequency Alfven waves appear to be the best candidates for coronal 
heating in open magnetic structures, such as coronal hole. The frequency of these waves 
determines the microphysics of the heating process. The low frequency (10“ 4 — 10“* Hz) 
Alfven waves can accelerate the solar wind by momentum transfer (first proposed by 
Alazraki k Courtier 1971. and Belcher 1971). and can heat the corona by viscous and 
resistive dissipation (first discussed by Barnes 1969). High frequency Alfven waves (i.e., 
10- 10 4 Hz range) will lead to the ion-cyclotron resonance which can preferentially heat the 
minor ions (e.g.. Dusenbery k Hollweg 1981. McKenzie k Marsch 1982), and was proposed 
as the heating and the acceleration of t he solar wind (e.g.. Axford & McKenzie 1992). 
However, it is not clear whether there is sufficient energy in the high frequency Alfven waves 
to account for heating of protons, and whether these waves can carry the energy from the 
magnetic network sufficiently high into the corona (e.g.. Cranmer. Field. k Kohl 1999). 

Recently, using the two-fluid 2.5D model we have shown that unresolved wave motion 
may account for the enhanced proton temperatures observed by UVCS (Ofman k Davila 
1999. Davila k Ofman 1999). However, the self-consistent effect of waves on the minor ions 
has not been considered. In the present study we use three-fluid electron-proton-ion plasma 
with a broad band Alfven waves source in the inillihertz frequency range and include the 
effects of Coulomb friction, important for minor ions. We consider the self-consistent effect 
of waves on the observed line width temperatures of protons and 0 5+ ions in a structured 
coronal hole. The results of this study are relevant for the future Solar Probe mission, that 
will measure the plasma properties near the Sun (Sittler et al 1997). 

UVCS observations indicate that 0 I,+ ions and other minor ions exceed temperatures 
of 10' K in corona] holes, and in streamers (Kohl et al. 1997). The UVCS observations 
show a striking difference in the appearance of coronal streamers in minor ion emission lines 
compared to hydrogen UY emission. These observations show that the emission from 0 0+ 
ions is stronger by an order of magnitude at the legs of streamers than in the central core 
of streamers (Figure 1). However, the brightness of the Ly-a emission peaks in the central 
core of streamers (Kohl et al. 1997, Raymond et al. 1997a). 

Several possibilities which can explain the observed difference in the 0 5+ emission and 
Ly-a in streamers are discussed by Noci et al. (1997). One of Noci et al. (1997) suggestions 
is that variation of minor ion abundance (i.e.. ion density relative to proton density) may 
occur due to Coulomb friction with outflowing protons (Geiss. Hirt. k Entwyler 1970) in a 
double-streamer geometry. Noci et al. (1997) concludes that only an abundance variation 
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of (.)°~ across the streamer is consistent with the observations of other minor ion (‘missions 
in streamers, such as Mg y ~ (625 A ) . Si fi_u (499A). and and Fe 1I+ (1242 A). Raymond et al. 
(1997a. b) suggested that gravitational settling in the core of streamers is the dominant 
mechanism for the depletion of 0 0+ there. However, the ID single-fluid time-dependent 
model of the settling and diffusion of ions in the streamer core predicts that the abundance 
in the streamer legs is less than in the core of the streamer, in conflict with observations 
(Raymond et al. 1997b). 

Early computational models of 0 ,J ^ ions in the solar wind did not include the full 
ion dynamics and treated the O'^ ions as test particles in one-dimensional solar wind 
flow (Geiss, Hirt, & Entwylcr 1970. Joselyn & Holzer 1978, Biirgi & Geiss 1986. Esser & 
Leer 1990). More recent multi-fluid studies of the flow of 0 5+ ions in the solar wind were 
confined to one spatial dimension (Li et al 1997. Raymond et al. 1997b). Although these 
models made important contribution to our understanding of the minor-ion dynamics in 
the solar wind, the effects of coronal magnetic and density structure due to streamers were 
not included self-consistentlv. 

Ion-cvclotron waves were suggested as the source of energy for the heating and the 
acceleration of the solar wind. There are various mechanisms that can generate ion-cyclotron 
waves in a plasma (e.g. ion temperature anisotropy, ion beams, etc.). However the dominant 
mechanism and its source for the generation of these high frequency waves in the solar 
corona is not well understood. If these waves are present, their direct signature could be 
detected in the future Solar Probe mission, that will measure in-situ plasma properties close 
to the Sun (Sittler et al 1997). 

It has been suggested that the observed 0 5+ temperature anisotropy is the result 
of ion-cyclotron heating by the appropriate power-law spectrum of ion-cvclotron waves 
generated lower in the corona (Cranmer, Field. & Kohl 1999). Cranmer et al. (1999) 
reports that the most probable perpendicular O 5-1 " speed is more than an order of magnitude 
larger then the most probable parallel speed at 3.5/iA. Due to the long integration time 
(several hours) of the UVCS observations and the line-of-sight contribution of several solar 
radii to the observed emission the observed velocity distributions are likely to contain both 
microscopic and macroscopic motions. The macroscopic motions may be due to unresolved 
magnetohydrodynamic (MHD) turbulence or low frequency MHD wave motions (Ofman & 
Davila 1997a). 



2. Three-Fluid 2.5D MHD Model Equations 


To investigate the dynamics of protons and ions in the solar wind we solve the 
three-fluid MHD equations with full ion dynamics (Braginskii 1965) in spherical geometry. 
The 2.5D equations are obtained by taking d/do — 0 in the 3D three-fluid MHD equations, 
i.e., by assuming azimuthal symmetry, and keeping all three components of the species 
velocities and fields. 

In the model we neglect electron inertia (since rn e rn p . where rn e is the electron mass 
and m p is the ion mass), relativistic effects (V < c), and assume quasi-neutrality of the 
plasma ( n e — n p + Zn r , where n r is the electron density, n p is the proton density, a* is the 
ion density, and Z is the charge state of the ions). We neglect the ion-cyclotron terms since 
the relevant time scales are orders of magnitude larger the the ion gyro-period. 

The normalized three-fluid MHD equations are 
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where the index k \ — p, i. (in Equation (5) k = e.p. /), and p is protons, i is ions, and e is 
electrons. In the present model we assume that the corona is nearly isothermal (i.e.. 1), 

do not include heat sources explicitly, and neglect the thermal forces in the momentum 
equations. The viscous forces assumed to be negligible. We also neglected a temperature 
equilibration term in the energy equation due to heat transfer between the fluids, since 
observations show that the 0 0+ temperature remains higher then the proton and electron 
temperature, and the protons will not be heated significantly by the 0 5+ due to the small 
abundance of these ions. In the above equations the Coulomb friction terms F^ cou / can be 
found in Braginskii (1965). Note, that qualitatively the Coulomb friction increases with the 
velocity difference between the species, and with the density. The friction decreases when 
the temperature is increased. 

We have used the following normalization: r — > r where R . is the solar 
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radius: t — » t/r A : \ ' — * \’/\' A \ B — ► B/ B {) : n imVtP -> /q.e.pAMfi 5 — w*/~ t the Lundquist 
number: 7 r = 4 -R'i/i/c 2 the resistive time scale, where v is the resistivity, and r is the 
speed of light: 7,1 = R. jV A the Alfven Lime scale: \’ A is the single-fluid Alfven speed: 

E Pp — ( k'bTfi'P'O I JJip ) j\ 2 the electron or proton fuller number: E, = (k h T Ki j/rn f )/\’ 2 \ the 
ion Euler number: F r ~ l f /? : / (G M . ) the Fronde number, where G is the universal 
gravitational constant and Al . is the solar mass: and the normalization constants 
a k = n p0 /n e0 /Ak + m P /m k + AkU^/n^ b = cB 0 f (4~en e0 R? ; IL* ) , where A k is the atomic 
mass of each species, and Z k is the charge state. 

We have used the following parameters in the three-fluid model of streamers: the 
proton temperature T pQ = 1.6 x 10 6 K: the electron temperature T e o — 1.6 x 10 6 K: the O 5 ^ 
temperature T a *+o = 12 x 10 6 K: the electron density n p0 — 10 H cnr ;i : and the O"^ density 
rii 0 = 6.2 x ID 4 cm " 4 : the magnetic field Bq — 2 G. These parameters are consistent with 
recent UVCS observations of streamers (Kohl et al. 1997. Li et al. 1998). Due to numerical 
limitations the 0 0+ temperature that we used is larger than the typical (P+ temperature 
(~ 5 x 10 h ) estimated for a coronal streamer. For the above temperatures the proton scale 
height is 0.14/? and the 0°+ scale height is 0.039 R.. Lower O s+ temperature results 
in smaller ion scale height, and nearly zero O 5 * 1 " density in the core of a streamer that 
poses a numerical challenge beyond the scope of this study. However, since the Coulomb 
friction increases for lower O 54 * temperature the results reported below should hold for 
T 0 s+ 0 = 5 x 10 6 K. 

To initiate the modeling of the coronal streamer we have used the quadrupole field 
given by B = J3 0 r _4 [(l - 3 cos 2 0)e r - sin 2 0e*]. The initial proton density and outflow 
was uniform in 9 and given by Parker's (Parker 1963) isothermal solution in r. The initial 
ion density was taken to be gravitationally stratified (with the gravitational acceleration 
corrected for the self consistent parallel electric field Raymond et ah 1997b), and uniform 
in 9. The initial uniform density at r - 1 R 0 was chosen to isolate the latitudinal density 
variations generated by the magnetic configuration above the limb as the simulation evolves 
to the steady state. The arbitrary initial state is chosen in order to reach the steady state 
solution in a reasonable computational time. The boundary conditions were imposed at 
r = with V> ?r > 0 (k = p, h), and with incoming characteristics approximated by zero 
order extrapolation. At r — r max: we have allowed for an outflow, and at 9 = 7r/2,7r we 
imposed = B 0 = 0, consistent with the quadrupole field geometry. The emphasis of 
this study is the slow solar wind, therefore, we did not include any additional sources of 
momentum or heating (such as waves) to produce the fast solar wind. 
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3. Slow solar Wind in Coronal Streamers: Numerical Results 

We solved the above equations explicitly with the fourth-order Runge-Kutta type 
method in time, and fourth order differencing in space with 160x500 grid cells (7)0 — 0.01. 
Sr = 0.01). A fourth order smoothing term was applied for stability. Two-fluid versions of 
this numerical approach have been used successfully for the modeling of the solar wind flow 
in coronal holes (Ofman k Davila 1999). 

The results of the calculations with the three-fluid code are shown in Figure 2 for the 
LVCS field of view. The model was run until a streamer with the appropriate sheet current 
was formed and reached a steady state. The current sheet has a finite thickness and is 
resolved in the code (the Lundquist number was set to 10°). The streamer proton density 
structure and magnetic field structure is in agreement with the single fluid streamer models 
(Pneuman. A: Kopp 1971. Steinolfson, Suess. k Wu 1982. Washimi. Yoshino. k Ogino 
1987, Cuperman. Ofman. k Dryer. 1990. Wang et al. 1993). The spatial variations of the 
proton density 6n pi and the 0 5+ density ()7q of a model coronal streamer at t = 16 hours 
are shown in Figure 2. 

It is evident that the 0 0i " density is enhanced in the open field line regions of the 
streamer compared to the core, and the proton density is larger in the core of the streamer. 
The density enhancement at the dose field boundary of the streamer is due to the 
compression by the Lorenz force due to the current sheet, calculated self-consistently in the 
model (the current sheet and the corresponding density enhancement is below resolution of 
present coronal observations). The proton density reached a steady state. The qualitative 
agreement between the model calculations of the proton and ion density structure, and the 
UVCS observations in Figure 1 is evident. Note, that in Figure 2 we model a streamer 
belt, while in Figure 1 the streamer is probably more localized in the longitudinal direction. 
Thus, 3D structure combined with line-of-sight integration of the observed 0 5+ emission of 
the "hollow" streamer may produce more localized legs. 

In the top panel of Figure 3 we show the magnetic field line vectors of the model 
streamer. The length of the streamline vectors represent the intensity of the magnetic 
field. The field lines are closed in the core of the streamer. The self consistent streamer 
field with current sheets results from the time-dependent 2.5D three-fluid calculation. The 
velocity streamlines are shown in the bottom panel. In the core of the streamer there is 
small downflow of O 0+ ions 20 km s -1 ), and no significant proton flow. Outside the core 
the flow is along the open field lines. The radial distance of the streamer cusp is determined 
by the magnetic field strength, geometry, and the plasma pressure (as in the single fluid 
models). 
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Fig. 1. Images of the UVCS 12 October 1996 observation of an equatorial streamer above 
the west limb in 0 0+ (top panel) and in Lyman-o (bottom panel) EUV emissions. The 
spatial scales are solar coordinates in arc min where the origin is at the center of the disk 
(Courtesy of SOHO/UVCS consortium. SOHO is a project of international cooperation 
between ESA and NASA). 

In the top panel of Figure 4 we show the 0-dependence of the radial velocity of ions 
(solid lines) and protons (dashed lines) at 1 . 8i? 0 (thick lines) and 5 R e (thin lines). Outside 
the streamer core it is evident that the outflow velocity of O r>+ is close to the outflow 
velocity of the protons. In the model we find that the outflow velocity is higher in the right 
leg of the streamer than in the left leg due to the asymmetry of the initial quadrupole field 
and the corresponding boundary conditions on the field at r = 1 i? 0 . The magnitude of the 
outflow speed in the open field regions is in agreement with the outflow velocity for the 
slow solar wind deduced from UVCS observations using the Doppler dimming effect at this 
height (Strachan et al. 2000). 

In the bottom panel of Figure 4 we show the 0-dependence of the proton and 0 5+ 
density. The 0 5+ abundance at the legs shows an order of magnitude enhancement 
compared to the core, in good agreement with the observations reported by Raymond et al. 
(1997a). The density spikes at the core boundaries are due to pressure balance with the 




Fig. 2. The spatial variations of the 0° T (top panel) and proton (bottom panel) densities 
in the model coronal streamer at t — 16 hours. 

current-sheet. They are not observed by UYCS, since the thickness of the current sheet is 
below resolution. The correlation between the outflow velocity, and the 0 5+ density is due 
to the effect of the Coulomb friction. The main differences in the detail between the model 
and observations stem from the difference in magnetic topology of the model streamer and 
the observed coronal streamer, and the possible latitudinal variation of the coronal density 
and temperature at the coronal boundary. 


4. Low-frequency MHD waves in solar coronal holes: Initial conditions, 

boundary conditions, and driver 

To model the dynamics of low-frequency MHD waves in the three-component plasma 
we use the following typical parameters: j p = 1.05, j e = 1.05, y 0 s + = 1.0. T p0 = 2 x 10 6 K, 
7 o 5 +o = 12 x 10 6 K. r„ 0 = 0.8 x 10 6 K, n e0 = 10 8 cm" 3 , n 0 s+ 0 = 8 x 10 4 cm' 3 . 5 = 10 4 . 
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Fig. 3. The magnetic field scaled by the square root of the total density (top panel) and 
the O 0+ velocity field (bottom panel) of the model coronal streamer at the time shown in 
Figure 1. The longest velocity vectors correspond to 144 km s -1 flows, and the longest 
magnetic filed vector correspond to 2.6 G. 

Dq = 7 G which gives l 4 = 1521 km s~b The abundance of 0 5+ relative to protons in the 
corona is about 8 x 10~ f) (Wimmer-Schweingruber et al. 1998). Due to numerical limitations 
we are using a higher abundance of 8 x 10~ 4 that is consistent with the photospheric 
abundance of oxygen relative to hydrogen (Anders & Grevesse 1989). However, this has 
no significant effect on the results, since the dynamics of electrons and protons is not 
significantly affected by the minor ions. The initial densities of the protons, ions, and 
electrons are (see Figure 5) 

rtp«,(0, r) = [l - (1 - n r ) e -[«W*)/*c] 4 ] 
n t o{0,r) = n p o(# ; l)n;i D (r), 

— n P o(ff. r) -f- 


n plD {r)/n r , 


( 6 ) 

(7) 

( 8 ) 
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Fig. 4. The ^-dependence of the velocity and the density at r — 1.8 R : (thick lines) and 
r ~ 5 R } (thin lines) for the model in Figure 1. The densities are scaled by the values at 
9 — 7T- Solid lines are for 0°A and the dashed lines are for protons. 

where n r = 0.3. the width of the low density region 9 C = 0.2. The ID densities n p iD(r), 
n i\D and velocities used to initialize the 2.5D calculations are obtained from tin? polvtropic 
steady state solution of the ID version of the code. 

The boundary conditions at r = 1 R# are incoming characteristics approximated bv 
zero order extrapolation. 

B r (l.9) = B r fr V ky0 = B $ = 0. T v = Tp 0 - T e = T e0 , T z = T l0 . (9) 

At r — f'max j and at 9 — 9 max we use open boundary conditions, with symmetry boundary 
conditions imposed at 9 = tt/2. 

The broad band Alfven wave driver is modeled by 

B4{t,8 J r = l) = -V d /V A . r F(L0) 

N 

F(t . 9) ~ ^ a, sin(^^ 4- r,(0)) 

i=i 


( 10 ) 

( 11 ) 
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Fig. 5. — The initial density profile of the protons, electrons, and O s+ ions in the model 
coronal hole. Note, that the electron and proton densities are nearly identical. The densities 
are scaled by 10 8 cm— 3. 

where a t — i 1/2 , = u>\ + {i — l)Au;, Ac u — — us\)/(N — 1). and 1^(0) is the random 

phase, and N is the number of modes. The frequencies in the broad band driver are in the 
millihertz range [u> 0 ,a;,v]. In the present study we use V d = 0 . 0151 % u; 0 = 0.013 Rad s -1 . 
u;,v = 10 %) > and N — 100 . The coronal boundary is divided into eight regions that produce 
the Alfven waves with a different 1^(0). 


4.1. Simulated emission line profile 

The calibrated observed line profile is usually given in terms of radiative intensity as a 
function of wavelength. In order to convert the wavelength to velocity the Doppler's shift 
relation is used 

A A v 

A 7 = C’ 


( 12 ) 
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Fig. 6. The thermal (solid line) and simulated (dashes) line profile at 4/7 .. The integration 
time is 1.7 hours. 

where Aq is the wavelength of the peak intensity of the emission line, r is the line of sight 
velocity, and A A is the wavelength shift due to the Doppler effect. 

The simulated line intensity is obtained by integrating the thermal emission shifted 
by the nonthermal wave velocity w(t) for the duration r (typically several hours) of the 
simulated observation and over the line of sight s : 


/(f) = f [ p e - {v - w(t))2/v B dt ds. (13) 

where u$ is the kinetic thermal speed of the particles, and the emission is weighted by the 
density of the species for the resonantly scattered line. In the present study we assumed 
that the line of sight contribution comes from about ±10 degree wedge in the center of the 
coronal hole. The effective temperature is obtained by determining the v\/ e of the simulated 
line width broadened by the nonthermal motion. An example of a simulated emission line 
profile broadened by nonthermal wave motions at 4 R s integrated over 1.7 hours is shown 
in Figure 6. 
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5. Waves in solar coronal holes: Numerical Results 

In Figure 7 we show the spatial dependence of the solutions at t = 109r.i = 13.9 hours. 
The o-eornponent of the proton and ion velocity, the Ohmic heating associated with the 
broad band Alfven waves, and the radial velocity of the protons are shown in the 0 — r 
plane. It is evident that, the Alfvenic fluctuations propagate into the coronal hole in protons 
and ions. We found that the ratio of the proton wave amplitude and the ion wave amplitude 
is about 3.2 = A/Z. This is expected from the ratio of the Lorenz forces exerted on the 
protons and on the ions in the momentum equation. It is also evident that the regions of 
largest Ohmic heating are associated with the highest amplitude Alfven waves, with phase 
mixing effects at the boundary between the low and high density regions of the coronal hole. 

The radial component of the proton velocity exhibits compressional fluctuations due 
to the nonlinear effects of the Alfven wave pressure, in agreement with single fluid models 
(Ofman & Davila 1997b. Ofman Davila 1998). On average, largest radial velocity occurs 
in regions of highest Alfven wave amplitude. The radial velocity amplitude predicted bv 
the model at 3 R® is about half the velocity deduced from UYCS observations at this height 
(Antonucci 1999). However, as Antonucci (1999) points out ’’...this velocity diagnostic 
depends critically on the modeling of the coronal properties along that line (line of sight)...” 
Thus, the reason for this discrepancy may be due to the choice of the magnetic field 
topology (i.e.. divergence), Alfven wave spectrum that affects the acceleration profile, and 
the other coronal parameters. Single fluid studies show that low frequency waves (<1 mHz) 
are needed to obtain the fast solar wind in coronal holes (Ofman & Davila 1998). However 
these low frequency waves accelerate the solar wind to fast wind speed beyond 10/? Q . 

The future Solar Probe mission will approach the sun to 4 i? 0 , and will be able to 
measure the variations in the particle velocities in-situ. To show' how these variations may 
look w r e plot in Figures 8-10 the temporal evolution of the V r and 1 / of protons (top panel) 
and 0 5+ ions (bottom panel). The location is r = 4 R P at 6 = 1.58 (solid line). 1.67 (short 
dashes), 1.76 (dashes), 1.85 (dot-dashes), 1.95 (3-dot-dashes). The velocities are scaled by 
the Alfven speed V A = 1521 km s -1 , and the time is scaled by the Alfven time, = 460 s. 

It is evident that the variations of V# of protons and ions are in phase. The amplitude 
of the Alfvenic fluctuations in ions is smaller than the proton amplitude by a factor of 
Z/A = 5/16. The temporal variations in V T of protons and ions (Figure 8) have global 
similarities. However, there are quantitative differences in the temporal variation due 
to different compressive effects in the proton and ion fluids. The driving source of the 
variations in V T are nonlinear compressions due to the low-frequency Alfven waves. 

The corresponding density variations of protons and ions at 4 i? ro are shown in Figure 9. 
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Fig. 7. I he spatial dependence of the ^-component of the proton and ion velocity, the 
Ohmic heating associated with the broad band Alfven waves, and the radial velocity of the 
protons are shown. The velocities are scaled by 0* = 1521 km s -1 . 


It is evident that the density variations occur in phase with the radial velocity variations, 
as expected from compressive fluctuations. The slow increase in the average ion density is 
due to the diffusion of 0 D+ ions to larger heliocentric distances due to Coulomb friction. 
Using the faster ID version of the code we have found that this process saturates after 
about 800r4. The time-averaged mass flux in the center of the model coronal hole at 4/2, ^ is 
2 x 10 g cm 2 s 1 , consistent with the mass flux at 1 AU determined from Ulysses data 
(Guhathakurta & Fisher 1998) for nearly radial expansion beyond 4 R Q . 

In Figure 11 we show the radial dependence of the effective proton and 0 5+ temperature 
in the central (low density) region of the coronal hole. The effective temperatures were 
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Fig. 8.- - The temporal evolution of the V r of protons (top panel) and O s+ ions (bottom 
panel). The location is r — 4 /? 0 at 9 = 1.58 (solid line), 1.67 (short dashes). 1.76 (dashes), 
1.85 (dot-dashes), 1.95 (3-dot-dashes). The velocities are scaled by V A = 1521 km s“\ and 
the time is scaled by the Alfven time, r A = 460 s. 
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Fig. 9. -- The temporal evolution of the proton density (top panel) and 0 0+ density (bottom 
panel). The location is r = 4 Rq at 0 = 1.58 (solid line), 1.67 (short dashes). 1.76 (dashes), 
1.85 (dot-dashes), 1.95 (3-dot-dashes). The densities are scaled by n e0 = 10 8 cm s“'\ and 
the time is scaled by the Alfven time, r A = 460 s. 
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Fig. 10. The temporal evolution of the V# of protons (top panel) and 0 5+ ions (bottom 
panel). The location is r = 4 R e at 9 = 1.58 (solid line), 1.67 (short dashes), 1.76 (dashes), 
1.85 (dot-dashes), 1.95 (3-dot-dashes). The velocities are scaled by V A = 1521 km s _1 , and 
the time is scaled by the Alfven time, r A = 460 s. 




obtained by integrating Equation (13) lor ])rot.ons and ions with tin 1 corresponding 
< Y,U) >. F he velocities were averaged in time and in space (0) to simulate the effects 
of 1 . i hour exposure of tlie LACS observations. Tin 1 dashed lines show the time-averaged 
kinetic temperature for the protons and ions that was obtained from the solution of the 
polytropic energy equation. The dotted line shows the electron temperature. 
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f ig. 11. The effective temperatures (solid lines) due to unresolved wave motions for protons 
and O ions obtained from the three-fluid model. The dashed lines are the time average 
temperatures obtained from the proton, ion, and electron poly tropic energy equations in the 
model. 

It is evident that the radial dependence of the effective proton temperature is in 
good qualitative agreement with UVCS observations (Kohl et al. 1997. Antonucci 1999). 
Due to the anisotropic velocity of the Alfven waves the apparent peak proton anisotropy 
is T e /f/Tp ~ 3. However, the Alfven waves can not explain the large perpendicular 
temperature of 0°+ ions reported by UVCS (Kohl et al. 1997, Antonucci 1999). The 
effective 0 5+ temperature in the model is only about 50% larger than the kinetic 
temperature, and an order of magnitude smaller than the perpendicular temperature 
deduced from UVCS observations. As we have seen above the amplitude of the heavy ion 
response to the Alfvenic fluctuation is reduced by the smaller Lorenz force compared to the 
force on the protons. Coulomb friction does not affect significantly the amplitude of the 
perpendicular ion fluctuations. Due to the small collision frequency with the protons the 
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main effect of the Coulomb friction is to increase the average ion out How speed and density 
on time-scales longer than the Alfven wave period (Ofman 2000). 


6 . Constraints on the O' Anisotropy due to the Ion-Cyclotron Instability 

In the present study we use linear kinetic theory and ID hybrid simulations to 
constrain the possible ion anisotropy for typical coronal hole conditions close to the 
Sun. We investigate the relaxation of the 0 1 ' ion temperature anisotropy for a range 
of ion parallel temperatures expected in a coronal hole. A 2-D hybrid code study of the 
anisotropic temperature oxygen ion instability with applications to the outer solar wind is 
in preparation (Gary et, al. 2000). 

We use the standard linear Vlasov theory for fluctuations and instabilities in a 
collisionless, homogeneous, magnetized plasma (Gary 1993). The protons and electrons 
are represented as Maxwellian velocity distributions, whereas the heavy ion species is a 
bi-Maxwellian with 7j_i/7|ji > 1. There are no average relative drifts between any two 
species, so the zero current condition is satisfied. Under these conditions, the fastest 
growing mode at sufficiently small % = Snn e k B T^/B^ is the electromagnetic heavy ion 
cyclotron anisotropy instability (hereafter the "heavy ion cyclotron instability"), where n, is 
the electron density, Bq is the background magnetic field, and k b is Boltzmann's constant. 
The maximum growth rate is at k x B 0 = 0; for this direction of propagation the linear 
properties of this mode are given by the dispersion equation stated in Gary and Lee (1994). 
We here present results from exact numerical solutions to this dispersion equation with no 
analytic approximations made for any species contribution. In each case we find that, on 
the given domain of (3 \\ <, the instability threshold condition can be fit to a curve of the form 

7]|i [(m p /m z )d |ij] ( ‘D ^ 14 

The fitting parameter Si changes with various parameters but is approximately of order 
unity, and a t ~ 0.40. 

Figure 12 shows the threshold anisotropies for the oxygen ion cyclotron anisotropy 
instability as functions of the oxygen ion scaling parameter (m p /mj)% for four different 
values of the normalized maximum growth rate y m /f l v . where Q p is the proton gvrofrequency. 
A stronger anisotropy corresponds to a faster growth rate, so although a, ~ 0.4 in each 
case, S t increases as 7 m/^p increases, just as S p increases with maximum growth rate at the 
threshold of the proton cyclotron anisotropy instability (Gary & Lee 1994). 

Other scalings of the linear instability threshold are not shown here, but may be 
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Fig. 12. Linear Vlasov theory results. The oxygen ion temperature anisotropy at four 
different values of the maximum growth rate of the electromagnetic oxygen ion cyclotron 
anisotropy instability as a functions of the scaling parameter (m p /m t ) Jj. Here the oxygen 
ions. 0 O+ . have T\\i/T\\ p = 16.0. ni t /m p — 16, and rq/n e ~ 0.0006. The fitting parameters of 
Equation (14) are: at y m /Q p = 0.001. 5 ? - - 1.46 and tv, - 0.42: at y m /Q p = 0.0025, 5 ?: = 
2.50 and = 0.38; at y m /Q p = 0.0050. S\ = 4.23 and o, : - 0.35: and at *f m /n p = 0.010, 
Si = 7.90 and a* = 0.34. 

summarized (under the assumption of 7 m /O p — constant) as follows: 5, generally decreases 
with increasing rii/n e and with increasing ionization state of the heavy ion (Gary et al. 
2000). In contrast, the threshold does not change significantly as T^/T p is varied. 

We are expanding our 2.5D MHD model into a full 3D description, which will allow to 
use more realistic model of a coronal hole and solar boundary conditions. For this purpose 
we are adapting our single fluid 3D MHD code and investigate various boundary conditions 
on the dynamics of waves in coronal magnetic structures (Ofman k Davila(2000)). Initially, 
we have developed the code on our SGI workstations. Next, the code has been parallelized 
by using the Message Passing Interface (MPI) to run on massively parallel computers, such 
as the T3E and on the Origin 2000. In the near future we plan to expand the code to model 
two- and three-fluid coronal plasma dynamics. 



7. Hybrid Numerical Simulations of the 3-component Plasma 


Wo use a one-dimensional hybrid code (Harned 1982. Ter as aw a et al. 1986. Winske 
& Omidi 1993. Agirn. \ inas. A' Goldstein 1995) to model a, collisionless. homogeneous, 
magnetized plasma. In the hybrid model the ions are represented as particles, while the 
electrons are described as a massless fluid to maintain quasineutrality of the plasma. This 
method allows one to resolve the ion dynamics and to integrate the equations over many 
ion-cyclotron periods, while neglecting the small time- and spatial- scales of the electron 
kinetic motions. The particle and field equations are integrated in time using a rational 
Runge-Kutta method whereas the spatial derivatives are calculated by pseudospectral FFT 
method. 

The electrons are represented as a massless fluid having zero average How speed. There 
are no initial relative drifts among the species. We use a periodic system in s with the 
background magnetic field Bq = xBq. We have used 256 grid cells to solve for the fields, 
with 100 particles per cell for each component and with grid size dx — 1.5 in units of proton 
inertial length, cf u; pp , where c is the speed of light, and uj vp is the proton plasma frequency. 
The system dimension L approximately equals four times the wavelength of the fastest 
growing mode. The integration time step was Q p A t = 0.05. 

We use the empirical values of the coronal parameters at 3.5/2.) reported bv Cranmer 
et al. (1999): n e — 3.55 x 10' 4 cm T T e — 7.2 x 10 5 K. 7]j p = 1.4 x 10 6 K. We assume that 
5 o(1/?q) - 5 G, and get Bq(3.5Rt) = 0.0655 G and — 6.3 x 10 2 Rad s' 1 . In this study 
we also assume that initially 7T p =7y p . Thus, we initialize the protons as a single isotropic 
Maxwellian velocity distribution, whereas the heavy ion species are initiallv a bi-Maxwellian 
with T ±i /T\\i = 50. 

The abundance of 0 5+ relative to protons in the corona at 3.5 R<-.) is about 8 x 10 -6 
(Wimmer-Schweingruber et al. 1998). Due to numerical limitations we are using higher 
abundance of 6 x 10” 4 . We have performed test runs with the coronal 0 5H " abundance, 
which requires much higher resolution, larger number of particles, and longer computation 
times, and found qualitatively similar results, with somewhat longer relaxation time and 
higher final anisotropy. To reduce the abundance to coronal values will require numerical 
effort beyond the scope of this study. 

We initiated the hybrid simulation with X[|j = Tj| p — 1.4 x 10 6 K. which corresponds 
t° $||i = 0.0413 (we vary below). We have run the simulation until the 0 5_t ~ anisotropy 
reached a quasi-steady state with final — 0.143. The 3\\i has increased due to the transfer 
of kinetic energy from the perpendicular to parallel direction of the ions. In Figure 13 we 
show the anisotropy of 0 5+ as a function of time. It is evident that the relaxation occurs in 



t rr i ~ 900U/. Hie anisotropy exhibits two main stages of evolution: (1) rapid relaxation 
on a typical timescale t rpi \ (h) slow decrease to an asymptotic value. 
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Fig. 13. The temporal evolution of the () ;>+ ion temperature anisotropy for typical coronal 
hole conditions. The initial parallel temperature was set to T\ t — T.\ p = 1.4 x 10 6 K, which 
corresponds to d\i = 0.0413. The final = 0.143. 

The rapid stage is characterized by an increase of the magnetic fluctuations evident 
in Figure 14. The increase of SB/Bl 2 is on the same timescale t r€i as the relaxation of 
the anisotropv The linear stage of the instability corresponds to an exponential growth of 
the magnetic perturbation. The average magnetic perturbation grows exponentially during 
the first 900 proton cyclotron periods, following an onset period of 300 proton cyclotron 
periods. At later times the instability nonlinearlv saturates. The magnetic fluctuation 
includes contribution from all unstable modes, with the dominant contribution from the 
fastest growing mode. 

In Figure 15 we show the velocity distribution of 0 0+ in the simulation at the end of 
the run in Figure 13. The distribution of (long dashes). V yA (solid line), and V sA (short 
dashes) is shown. The (lotted line is the corresponding Maxwellian. The V x%i is parallel, and 
the \ ya and l zA are perpendicular to the background magnetic field. It is evident that the 
anisotropic distribution of the 0 o+ velocity is close to bi-Maxwellian with the perpendicular 
component broader than the parallel in agreement with the anisotropy at the end of the 
run in Figure 13. The fluctuations in the distributions due to the statistical noise do not 
not significantly affect the results. 

The relaxation of the anisotropy and the increase of the magnetic fluctuation occur 
due to the ion-c\clotron instability that develops in the The rapid phase corresponds 
to the exponential growth phase of the instability, while during the asymptotic relaxation 
stage the instability has nonlinearlv saturated. Linear Vlasov theory, and nonlinear particle 
and hybrid simulations show that in collisionless plasma the relaxed velocity distribution 
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Fig. 14.— The temporal evolution of the perturbed relative magnetic field fluctuation 
|dB/B| 2 generated by the ion-cvclotron instability of the 0 5 ~ ions. To remove statistical 
high-frequency noise a running average smoothing was applied to this quantity. Note, 
that the saturation time of the fluctuations corresponds to the relaxation time of the ion 
anisotropy. 
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Fig. 15. The velocity distribution of O s+ in the simulation at the end of the run of 
Figure 13. The distribution of V' x . t (long dashes), V y4 (solid line), and (short dashes) are 
shown (N is the number of the particles). The dotted line is the corresponding Maxwellian. 
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does not return to a fully isotropic condition (Cuperman. ()fnian. N Dryer 1986. tiarv et 
al. 2000). 
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Fig. 16 — The results of the parametric study of the hybrid simulation. The scaling of the 
relaxed T± t /T\^ — 1 with the final (full circles). The scaling of the relaxation time. t re [ , 
with the initial % (circles). Both quantities scale as .T- 0 ' 11 . The *’x' J, s mark the values 
- 1 at t— 0. 

In Figure 16 we show the dependence of t re i and the relaxed anisotropy on In this 
study we have varied the initial in the range 0.041-0.577. This corresponds to T\n in the 
range (1.4 - 20.0) x 10 6 K. We found that t re i and the final T u /T\ {i - 1 scale as J^ Al in 
agreement with the linear theory (Equation 14). Note, that we plot the scaling of t re i with 
the initial 3\\%^ while the scaling of 7j_i/Tj| z - — 1 is with the final 3\^. The ,l x J ”s denote the 
anisotropy at t— 0 for each case. We find that decreasing the initial Tii, leads to larger final 
anisotropy of the 0°"F and to longer relaxation times that correspond to smaller growth 
rate of the ion-cvclotron instability. 




8. Analysis 


We developed the three-fluid 2.5D MHD numerical model to investigate the acceleration 
of the slow and fast solar wind, and the effect of waves on protons and minor ions. According 
to this model the main mechanism responsible for the observed structure of streamers 
in 0° + is the Coulomb friction between the () 5 ~ ions and the outflowing protons. The 
friction leads to the increase in the density and the outflow speed of the minor ions in 
the open field regions. The same qualitative enhancement is expected for Fe 11 * and 
other minor ions in these regions. In order to verify that the Coulomb friction is the 
physical mechanism that, leads to the enhanced 0° + abundance in streamer legs we have 
repeated the above calculation without the Coulomb friction terms, and found much smaller 
abundance enhancement of 0 0-h (order of few percent) in the legs of streamers, and an 
order of magnitude smaller 0 o+ outflow speed in the open field regions. Other possible 
mechanism that may contribute to the enhanced abundance in the legs of streamers 
are preferential heating of the ions in these regions, and variations in the ion abundance at 
the coronal boundary of the legs. These possibilities are left for future studies. 

Based on these results, and on recent SOHO/UVCS observations of the solar wind in 
solar coronal streamers we conclude that the observed enhanced minor ion emission in the 
legs of streamers is the signature of the slow solar wind outflow regions. The identification 
of these regions allows to trace the slow solar wind sources low above the limb, and into 
the extended corona. This will help to understand the origins and the composition of the 
slow solar wind, and the effect on the interplanetary environment. Using the three-fluid 
2.5D MHD model we study the role of broad-band low frequency MHD waves and coronal 
structures on the EUV line profiles. We calculate the apparent proton and ion temperatures 
from the simulated EUV emission lines. We find that for reasonable coronal parameters 
and Alfven wave spectrum (i.e., consistent with nonthermal line broadening near the limb) 
the wave motions in the line of sight can account for the observed Ly-o: line broadening, 
but can not explain the broad OVI lines. 

We find with our three-fluid model that the UVCS observations of broad Ly-u lines 
suggest the presence of low-frequency Alfven wave in coronal holes that are responsible 
for the enhanced proton temperature and anisotropy. Thus, these waves accelerate the 
protons in coronal holes in agreement with our earlier single fluid models (Ofman k Davila 
199 /a, Ofman k Davila 1997b, Ofman k Davila 1998). Our three-fluid model suggest 
that the observed increase of the line-width proton temperature is due to temporally and 
spatially unresolved low-frequency wave motions of the protons. Thus, the MHD waves may 
contain a significant fraction of the solar wind energy in the lower corona, and therefore, 
may play an important role in the acceleration and heating of the fast solar wind. 
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However. tin* low frequency Alfven wave spectrum (an not product 1 the broad ()'~ lint's 
observed by the l \ CS. According to tin 1 three-fluid model the perpendicular ()*~ velocities 
art 1 about factor of 5 smaller than the velocities reported by Kohl et al. (1997). Due to the 
low abundance, low collision rate, and low charge to mass ratio of the 0 5+ ions compared 
to protons it is possible that high frequency waves play a more important role in heating 
these ions than in proton heating. Thus, resonant absorption of ion-cyclotron waves, may 
be responsible for the high minor ion temperature and anisotropy deduced from L : YCS 
observations in agreement with other studies (e.g.. Cranmer. Field, k Kohl 1999. Ofman. 
Vinas. k Gary 2001). 

We find that the large 0°^ ion anisotropy deduced from U\CS observations is unstable 
to an electromagnetic ion-cyclotron instability. We use the scaling law for instability 
threshold obtained from linear \ lasov theory to relate the heavy ion temperature anisotropy 
to the parallel heavy ion beta. The linear theory relies on the assumption that the ion 
distribution function is bi-Maxweilian. which is reasonable for the inner corona. Obviously, 
the linear theory does not account for the self-consistent interaction between the waves and 
the particles, and the saturation of the instability. 

The hybrid simulation allows the self-consistent evolution of the velocity distribution of 
the ions (e.g.. bi-Maxwellian distribution is not imposed at t > 0). It includes the nonlinear 
self-consistent description of the wave-particle interaction, and the nonlinear saturation of 
the instability. The simulations show that these nonlinear effects constrain the values of the 
ion anisotropy for a range of coronal parameters. Furthermore, the simulated constraint is 
in agreement with the instability threshold scaling derived from linear theory. 

For typical coronal hole conditions at 3.57? 0 and isotropic protons, an initial 0 5+ 
anisotropy of 50 decreases by an order of magnitude within ~300-900 proton cvclotron 
periods (3-9 s for coronal values of B). The relaxation time and the final anisotropy 
decreases when the initial T.\i is increased. Thus, in order to sustain the anisotropy of 50 
the 0 5+ ions must be heated and the resonating ion-cyclotron wave spectrum must be 
produced on these timescales. The results of this study are applicable to the future Solar 
Probe in-situ velocity distribution observations close to the Sun. 
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